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Using the large hadroproduced charm sample collected in experiment 
E791 at Fermilab, we have measured ratios of branching fractions 
for the two-body singly-Cabibbo-suppressed charged decays of the D°: 

t(d^ k4+) = a109 ± a003 ± a003 ' rffi -> = ao40 d= 

0.002 ± 0.003, and K 1 K -) = 2.75 ± 0.15 ± 0.16. We have 

1 [D — > TT^TT 

looked for differences in the decay rates of D° and D° to the CP eigen- 
states K + K~ and ir + ir~, and have measured the CP asymmetry pa- 
rameters A C p{K + K-) = -0.010 ± 0.049 ± 0.012 and A C p{it + it~) = 
-0.049 ± 0.078 ± 0.030, both consistent with zero. 

PACS: 11.30.Er; 13.25.Ft; 14.40.Lb 



The measured world average for the ratio of the branching fractions of the 
two-body singly-Cabibbo-suppressed (SCS) charged decays of the D° meson 

is ^7 „n 1 _, J = 2.86 ±0.28 pjl. Models including final state interactions 

1 [JJ — > 7T^7T ) LJ1 

0-0, penguin diagrams 0, QCD sum rules 0, and non-perturbative algebraic 
approaches [0 have been proposed to explain the experimental value observed 
for this ratio. Precise measurements of this ratio can help to differentiate 
among models, and can also aid in our understanding of the Standard Model 
predictions for D°-D° mixing via long-range mechanisms, which require SU(3) 
symmetry-breaking to be non-zero |8HiO . 



To date CP violation has been observed only in the neutral kaon system. In 
the Standard Model this violation is a consequence of a complex amplitude in 
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the Cabibbo-Kobayashi-Maskawa matrix. In this model, strange quarks couple 
to the top quark in diagrams with internal loops, leading to observable CP 
violation. The comparable diagrams in the charm sector have bottom quarks 
in the internal loops, resulting in very small Standard Model contributions 
to CP violation. Thus the charm sector may be uniquely sensitive to physics 
outside the Standard Model, at the 1CT 3 level [[ll|]. Currently, the measured 
limits for charged and neutral Z)'s are at the level of (5-10)% from experiments 
E687 p|,|ll and E791 [111 at Fermilab, and CLEO PLpl. 



In this paper, we first report measurements of the branching fractions for the 
SCS decays D° — > K + K~ and D° — > 7r + 7r~ relative to the Cabibbo-favored de- 
cay D° — > K~n + , based on a total sample of about 42,000 fully reconstructed 
two-body D° decays. Second, we present a search for CP violation in the 
SCS decay of the D° using a sample of 14,225 Z)°'s tagged through the decay 
D* + -> D°ix + . Throughout this paper, the CP-conjugate states are implicitly 
included unless otherwise noted. At the level of sensitivity of this experiment, 
any direct CP violation would indicate physics outside the Standard Model. 



The current results are based on data accumulated by experiment E791 in a 
500 GeV 7r~ beam during the 1991/92 Fermilab fixed-target run. E791 was 
the fourth in a series of charm experiments performed in the Fermilab Tagged 
Photon Laboratory. The E791 spectrometer ]17| was an open geometry detec- 
tor with 23 planes of silicon microstrip detectors (6 upstream and 17 down- 
stream of the target), 35 drift chamber planes, 10 proportional wire chambers 
(8 upstream and 2 downstream of the target), two magnets for momentum 
analysis, two large multicell threshold Cerenkov counters for charged particle 
identification, electromagnetic and hadronic calorimeters for electron/hadron 
separation as well as for online triggering, and a fast data acquisition system 
that allowed us to collect data at a rate of 30 Mbyte/s with a 50 /zs/event 
deadtime. The target consisted of a 0.6-mm thick platinum foil followed by 
four 1.5-mm diamond foils. Each target center was separated from the next 
by about 1.5 cm, allowing observation of charmed-particle decays in air with- 
out background from secondary interactions. The very open transverse-energy 
trigger was based on the energy deposited in the calorimeters and was highly 
efficient for charm events. Over 2 x 10 10 events were recorded during a six- 
month period. 



The two-body decay sample was selected based upon several criteria. A decay 
track candidate must have made a large contribution to the \ 2 °f the event 
production vertex fit when included in that fit. The significance of the mea- 
sured separation (in the beam direction) of the candidate decay vertex from 
the production vertex had to be > 8a, where a is the error on the measured 
separation of the two points. The momentum component of the D° candidate 
transverse to the line connecting the production and decay vertices had to be 
less than 0.40 GeV/c. The sum of pf of the decay tracks, with p t measured 
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relative to the direction of the D° candidate, was required to be greater than 
0.52 (GeV/c) 2 . Finally, the decay vertex had to be located well outside the 
target foils. 

Y(D° — * K + K~) 

To improve the statistical significance of the ratio J;, , _J , we re- 

1 (D U — > 7T + 7T ) 

quired the daughter tracks to have K or tc signatures in the threshold Cerenkov 
counters. However, the efficiency of our identification depends on particle mo- 
mentum. We have studied these particle identification efficiencies for the data 
sample D° — > K~n + as a function of momentum, transverse momentum with 
respect to the beam direction, charge, and particle type. Based on this study, 
we required the particle momenta to be in the range 6 to 80 GeV/c. We 
then corrected for the particle identification efficiency by weighting each D° 
candidate by the inverse of the product of the two particle identification ef- 
ficiencies. Using this procedure we have weighted each track according to its 
p and p t , rather than applying one global factor to the final results. We cal- 
culated the statistical errors in the weighted signals by scaling the statistical 
errors of the unweighted signals by the ratios of the weighted to unweighted 
sample sizes. Figure [1] shows the Cerenkov-weighted invariant mass distribu- 
tions for our data. The peaks to the right of the signal region in the K + K~ 
distribution and to the left of the signal region in the 7r + 7r~ distribution are 
due to misidentified Kn events, whereas the peak to the left of the signal re- 
gion in the K + K~ distribution is due to misidentified Knir events. A binned 
maximum-likelihood fit was performed for each distribution with the signal 
function assumed to be Gaussian. The backgrounds were fit to the following 
functions: K~tt + - a third-order polynomial (Figure H left); K + K~ - a linear 
term plus a Breit-Wigner function for misidentified Knn and a half-Gaussian- 
half-Breit-Wigner function for misidentified Ktx (Figure [l] middle); 7r + 7r~ - an 
exponential term plus a half-Breit-Wigner-half-Gaussian for misidentified Ktt 
(Figure |] right). Cross-hatched areas are the estimated backgrounds under 
the signals. 



We have used a Pythia-based Monte Carlo (MC) [18| which incorporates a 



full detector simulation to correct for detector acceptance effects other than 
Cerenkov identification. The efficiencies determined from the Monte Carlo sim- 
ulation and the number of signal events extracted from both the unweighted 
and weighted mass distributions are given in Table [I]. 

Table 1 

Sample sizes and average geometric acceptances for three decay modes of the D°. 





D° -> K~n + 




K~ 


D° -> 7T + 7T- 


Unweighted Sample 


36955 ±217 


3317 ± 


84 


2043 ± 95 


Weighted Sample 


63177 ±371 


6845 ± 172 


2521 ± 117 


Average Acceptance 


3.21% 


3.18% 




3.22% 
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1.8 2 2.2 1.8 2 2.2 1.8 2 2.2 



M(KV) (GeV/c 2 ) M(K + K") (GeV/c 2 ) M(ttV) (GeV/c 2 ) 

Fig. 1. Invariant mass of Cerenkov- weighted candidates. The solid lines correspond 
to a binned maximum-likelihood fit to a signal given by a Gaussian distribution plus 
different background functions given by: D° — > K~ir + - a third-order polynomial 
(left); D° — > K + K~ - a linear term plus a Breit-Wigner function for misidentified 
Knir and a half-Gaussian-half-Breit-Wigner function for misidentified Ktt (mid- 
dle); and D° — > tt + tt~ - an exponential term plus a half-Breit-Wigner-half-Gaussian 
function for misidentified Kit (right). Cross-hatched areas are our estimates of back- 
grounds below the signals. 

We calculate the ratios of branching fractions from the numbers in the last 
two rows of Table [I] and present them together with previous measurements in 
Table ||. The quoted systematic errors are the quadrature sums of systematic 
uncertainties from the following sources: relative efficiencies for the selection 
criteria, fitting functions, MC production models, and the Cerenkov particle 
identification weighting procedure. Table § presents the contribution of each 
of these to the total systematic uncertainty, expressed as a percentage of the 
statistical uncertainty. The "Selection Criteria" entry reflects uncertainties in 
the details of the MC modeling of the experimental acceptance. The "Fitting 
Functions" uncertainty corresponds to various estimates of the background, 
especially for tttc. We note that reflections do not contribute directly in 
the KK or ttit signal regions. The shapes of the Kir and Kinr reflections have 
been studied in detail, and the same fitting function yields good fits to both 



5 



MC and data. Uncertainties resulting from the MC production model were 
investigated by changing PYTHIA default parameters to agree with a study 
of D ± production in E791fll9l 



Table 2 

D° — » K + K~ and tx + ix~ relative branching ratio measurements compared with 
previous experiments. 



Year 


Group 


r(L> u -> K+K~) 

r(D u -> K-7T+) 


T(D [) -> tC 

r(L> u -> k 


-7T-) 

-7T+) 


T(D U -> 7T + 7T-) 


1997 


E791 


0.109 ±0.003 ±0.003 


0.040 ±0.002 


±0.003 


2.75 ±0.15 ±0.16 



1979 
1984 
1989 
1990 
1991 
1992 
1993 
1994 
1996 



Mark II [[21]] 
Mark III p2| 
ARGUS H 
CLEO P] 



E691 f24 
WA82 |25 
CLEO |B| 
E687 m 
PDG El 



0.113 ±0.030 0.033 ±0.015 3.4 ±1.8 

0.122 ±0.018 ±0.012 0.033 ±0.010 ±0.006 3.7 ±1.3 

0.10 ±0.02 ±0.01 0.040 ±0.007 ±0.006 2.5 ±0.7 

0.117 ±0.010 ±0.007 0.050 ±0.007 ±0.005 2.35 ± 0.37 ± 0.28 

0.107 ±0.010 ±0.009 0.055 ±0.008 ±0.005 1.95 ± 0.34 ± 0.22 

0.107 ±0.029 ±0.015 0.048 ±0.013 ±0.008 2.23 ± 0.81 ± 0.46 

0.0348 ± 0.0030 ± 0.0023 

0.109 ±0.007 ±0.009 0.043 ±0.007 ±0.003 2.53 ± 0.46 ± 0.19 

0.113 ±0.006 0.0396 ±0.0027 2.86 ± 0.28 



Table 3 

Contributions to the systematic uncertainty for each of the measured branching 
ratios. The systematic uncertainties are expressed as a percentage of the statistical 
uncertainty for the corresponding branching ratio 

r(£> u -> K + K 

TW 



K~tx a 



TX ' IX 



Tp 1 



K + K' 



K~tx- 



TX 1 TX 



Selection Criteria 


90% 


100% 


40% 


Fitting Functions 


60% 


70% 


90% 


MC Production Model 


30% 


30% 


20% 


Particle ID Weighting 


30% 


30% 


30% 


Total 


116% 


129% 


105% 



We now present results of a search for CP violation in the SCS decays of the 
D° (or D°) beginning with the same data sample as above. Assuming CP 
conservation in the strong decay of the D*, we identify the meson as either 
D° or D° by tagging it using the charge of the slow tx from D* + — > D°tx + and 
D*~ — > D°7x~ . We required the mass difference between D* and D to be in the 
range 143 - 148 MeV/c 2 and the distance of closest approach of the bachelor 
pion to the primary vertex had to be < 120 /im. Because D* tagging reduces 
background so strongly, we relaxed the requirement on the sum of pf of the 



6 



D° decay tracks relative to the direction of the parent D° from 0.52 (GeV/c) 2 
to 0.4 (GeV/c) 2 . For this analysis, we used the same Cerenkov identification 
criteria as for the branching ratio measurements, but we did not weight the 
events. 

The signal for CP violation is an absolute rate difference between decays of 
particle and antiparticle to charge-conjugate final states / and /, where / = 
K + K~ or 7r + 7r~: 

_ r(D->/)-r(p->7) m 

CP V(D -»•/) + Y(D -> /) {> 

In hadroproduction, D and D mesons are not produced equally. Therefore we 
normalized our signals to the Cabibbo-favored D° — > K~tt + and D° —> K + tt~ 
signals and measured 



N(D° 


-/) 


iV(D° 


-/) 


N(D [> -> 




-) iV(.D u 




N(D° 


-/) 


N(D° 


-7) 


N(D° -> 









(2) 



where, for each channel, iV = en is the observed number of events, n is the 
produced number of events, and e is the efficiency of our detector and analysis 
procedure. To the extent that 

e(D° - f) = <D° - J) ^ 



the measured quantity (Eq. 2) is A C p- We verified the validity of Eq. 3 for this 
analysis in two stages. First, we used real data to determine that the Cerenkov 
identification efficiencies for kaons and pions are independent of charge in any 
given momentum range from 6 to 80 GeV/c. Then, we used our Monte Carlo 
simulation of the experiment to determine that the geometric acceptances are 
independent of charge. 

Note that any CP asymmetry from interference between mixing and tree- 
level diagrams will not cancel through the D° — > K~n + normalization. Thus 
our normalized Aqp is not a direct CP asymmetry parameter, but rather a 
measure of combined direct and indirect CP asymmetries ||2U|| . An implicit 
assumption in this analysis is that there is no measurable CP violation in the 
Cabibbo-favored decays of the D°. 

Figure |2| presents mass plots for the candidate D° and D° decays to K-k, KK, 
and 7T7T, including our fits to the distributions. We used the same fixed central 
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1.7 1.8 1.9 2 2.1 2.2 



M(K"tt + ) (GeV/c 2 ) 




M(K + K") (GeV/c 2 ) 




1.7 1.8 1.9 2 2.1 2.2 

M(tt + tt") (GeV/c 2 ) 




1.7 1.8 1.9 2 2.1 2.2 



M(K + tt") (GeV/c 2 ) 











D° - 


> K + K" 




333 ± 21 


1 1 1 I 1 1 


ill i i 1 





1.7 1.8 1.9 2 2.1 2.2 



M(K + K") (GeV/c 2 ) 




M(tt + tt") (GeV/c 2 ) 



Fig. 2. Invariant mass of D° candidates (left column of figures) and D° (right col- 
umn of figures) tagged with the decay D* — > D°ir, for the decay D° —* i"C _ 7r + (top), 
D° —> K + K~ (middle), and D — > tt + tt~ (bottom). The numbers presented are 
the sample sizes as calculated from Gaussian fits to the signal regions, and the 
cross-hatched areas are our estimates of backgrounds below the signals. 



mass value and signal width for D° and D° fits. The distributions fit the data 
well, and the integrals of the Gaussian fits to the signals were used to calculate 
the observed numbers of decays. Our asymmetry results are listed in Table (| 
along with previous measurements. We see no evidence of CP violation. 
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Table 4 

The measured asymmetry parameters Acp for K + K~ and tt + it~ decay modes of 
the D° and their 90% confidence-level limits, compared with previous experiments. 



1 Cell 


Group 




T iiTni"!" n "1" Qflvr\ r'OTi "h H <^ti r*fi 1 p\rp 1 

±J 11111 L CL u U \J /(J ^VjllllLldlL-C 1C Vv/1 


1997 
1997 


E791 
E791 


-0.010 ±0.049 ±0.012 
-0.049 ±0.078 ±0.030 


-9.3% < A CP (K+K-) < 7.3% 
-18.6% < A CP (n + n-) < 8.8% 


1991 
1994 
1995 


E691 [ 
E687 [ 
CLEO I 


24 
13 


| 0.20 ±0.15 
| 0.024 ± 0.084 
26f 0.080 ±0.061 


A CP {K + K-) < 45% 
-11% < Acp(K+K-) < 16% 
-2.2% < A CP (K + K-) < 18% 



In summary, we have measured the following ratios of branching fractions for 

r(£>° -> AT+AT 



the charged two-body decays of the .D : 



0.003, 



r(z> 



AT-7T- 



r(D° -> K-n^ 
0.040 ± 0.002 ± 0.003, and 



^ = 0.109 ± 0.003 ± 
r(L>° -> AT+A'-) 

r(D° ^ tt+tt 



2.75 ± 0.15 ± 0.16. We also find the 90% confidence-level intervals on the 
decay asymmetries of D° and D° to the CP eigenstates K + K~ and it + it~ 
to be -9.3% < A CP (K+K~) < 7.3% and -18.6% < A CP (n+n-) < 8.8%. 
We find no evidence for CP violation in these decay modes. This is the first 
reported result of a search for CP violation in the decay mode D° — > 7r + 7r~. 

We gratefully acknowledge the assistance of the staffs of Fermilab and of all 
the participating institutions. This research was supported by the Brazilian 
Conselho Nacional de Desenvolvimento Cientifico e Tecnologico, CONACyT 
(Mexico), the Israeli Academy of Sciences and Humanities, the U.S. Depart- 
ment of Energy, the U.S. -Israel Binational Science Foundation, and the U.S. 
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